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We discuss several issues important for experimentally observing Efimov physics in ultracold 
quantum gases. By numerically solving the three-boson Schrodinger equation over a broad range of 
scattering lengths and energies, and by including model potentials with multiple bound states, we 
address the complications of relating experimental observations to available analytic expressions. 
These more realistic potentials introduce features that can mask the predicted Efimov physics at 
small scattering lengths. They also allow us to verify that positive and negative scattering lengths 
are universally connected only across a pole, not across a zero. Additionally, we show that the 
spacing between Efimov features for the relatively small scattering lengths accessible experimentally 
fail to precisely follow the geometric progression expected for Efimov physics. Finally, we emphasize 
the importance of the short-range three-body physics in determining the position of Efimov features 
and show that theoretically reproducing two-body physics is not generally sufficient to predict three- 
body properties quantitatively. 

PACS numbers: 34.50.-s,34.10.+x,21.45.+v,05.30. Jp 
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I. INTRODUCTION 

The importance of three-body collisions in determin- 
ing the lifetime and stability of ultracold atomic gases has 
been underscored in recent experiments Near a Fes- 
hbach resonance, for instance, the two-body s-wave scat- 
tering length a can be made much larger than the char- 
acteristic range r$ of the interatomic interactions and, 
due to the strong dependence of the ultracold three-body 
collision rates on a 0, the loss processes can dominate, 
burning the sample completely or, in other cases, lead- 
ing to unexpected stability when such collisions are sup- 
pressed [3J- In addition to these practical considerations, 
fundamental effects associated with Efimov physics also 
appear in the strongly interacting limit ja\ 3> tq. In 
the early 1970's Vitaly Efimov predicted [J], in the con- 
text of nuclear physics, the existence of a large number 
of exotic weakly bound trimers that appear whenever 
\a\ 3> tq. Only recently, though, has the first experi- 
mental evidence of such Efimov states been obtained in 
ultracold gases [f|, representing a landmark step for a 
new generation of experiments designed to explore the 
rich world of few-body physics. 

In ultracold quantum gases, Efimov physics is revealed 
through its influence on collisions rather than through 
observation of the states themselves. In this case, Efimov 
physics can be observed by measuring the time evolution 
of the atomic or molecular densities and extracting the 
collision rates [5]. Efimov physics appears as a series of 
minima or maxima in these rates as a function of a with 
the features separated by the multiplicative factor e*' s ° , 
where sq = 1.00624 for identical bosons but can vary 
greatly for other systems @ . Many of these features are 
associated with the creation of a new Efimov state in the 



system. 

In this paper we discuss some of the issues that must be 
considered when attempting to establish connections be- 
tween analytical predictions for Efimov physics (see, for 
instance, the reviews [9j, [lC§ and Ref. [2j) and experimen- 
tal data. These analytic expressions are derived as ex- 
pansions about T—0 and |a|/ro = oo, typically including 
only the leading term of the expansion. The limitations 
of these expressions have been explored both in temper- 
ature 0] and in scattering length @, [!]• In particular, 
it has been shown that the \a\ 3> ro limit is approached 
only very slowly for many key properties of three-body 
systems, leading to complications in applying the ana- 
lytic expressions to experiment since, unfortunately, ex- 
periments cannot easily access this limit. Appropriately 
corrected analytical formulas, although potentially uni- 
versal themselves, are not yet known [3, [13]. Perhaps 
surprisingly, another class of complications arises from 
the fact that some experiments — even though they are 
performed at seemingly "ultracold" temperatures — ac- 
tually lie outside the range of validity 0] of the zero tem- 
perature analytical expressions currently available. Fi- 
nally, there are still other complications that can arise 
for realistic interatomic potentials that have not been 
incorporated into the analytical expressions. These com- 
plications include, for instance, the finite range of the 
potentials, the existence of higher angular momentum 
two-body states, and the multichannel nature of the two- 
body interaction. 

While progress is being made to find more general an- 
alytic expressions pj, Oj, [H HH, [H, H3, 13] , direct nu- 
merical solution of the Schrodinger equation currently 
remains the most reliable approach for non-zero colli- 
sion energies, finite range potentials, and finite a. We 
present our numerical results for three-boson recombina- 
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tion, B+B+B—*B2+B, using the adiabatic hyperspher- 
ical representation 0, 0] for a wide range of collision 
energies. Our calculations show that if the \a\ ^> r con- 
dition for Efimov physics is not fully satisfied, then the 
spacing between Efimov features can be quite far from 
the predicted e^ s °. Moreover, we show that if the two- 
body model has several bound states, then additional 
resonant effects can appear, complicating the interpreta- 
tion of Efimov features. We also demonstrate that if a 
is changed from +oo to — oo by crossing a = 0, then the 
relation between Efimov states that occur for a > and 
a < is not universal, unlike the case when a changes 
instead through \a\ = oo across the resonance [loj]. To 
illustrate these points, we present numerical results for 
recombination using model interactions chosen to match 
the recent Innsbruck Cs experiment [5J. 

Finally, we discuss the importance of the short-range 
three-body physics in determining the precise position 
of the Efimov features. In our work, this physics ispa- 
rameterized by a short-range phase (also in Ref. 
whereas effective field theory treatments have various pa- 
rameterizations generally equivalent to a momentum cut- 
off [1,0. This three-body parameter cannot, however, be 
determined from knowledge of the near-threshold two- 
body physics. Consequently, fitting the near-threshold 
two-body observables [13, EH alone is not sufficient to 
accurately predict the position of the Efimov features, 
no matter how complete. To make this point clear, we 
show that including a non-pairwise-additive three-body 
interaction [22L |23(, known to be present in every tri- 
atomic system, changes the positions of the Efimov fea- 
tures, which is equivalent to changing the short-range 
three-body phase [jjj • So, even if the "exact" two-body 
potentials were used instead of just the near threshold fit, 
the positions of the Efimov features could not be quan- 
titatively predicted. 



the three-body problem can be described by the hyper- 
radius R and five hyperangles, denoted collectively by f2. 
The hyperangular part of the Schrodinger equation de- 
termines the relative motion of the three bodies, and the 
hyperradial part determines the overall size of the system. 
The five angular coordinates are chosen to be the Eulcr 
angles (a, /3, and 7), specifying the orientation of the 
plane defined by the three particles relative to the space- 
fixed frame, plus two hyperangles (8 and <p) defined as a 
modified version of Smith- Whitten coordinates [13, HH| ■ 
The Schrodinger equation in hyperspherical coordi- 
nates can be written in terms of the rescaled wave func- 
tion tp = i? 5 / 2 ^ as (in atomic units) 



1 d 2 

2jldR? 



H ad {R,n) 



r/>(R,Sl) = EtP(R,Q), (1) 



where \i is the three-body reduced mass (/1 = m/V3 for 
three identical bosons of atomic mass m) and E is the 
total energy. The adiabatic Hamiltonian H a d in Eq. ((T|) 
is given by 



H ad {R,fl) = 



A 2 + 15/4 
2^iR? 



V(R,0,<p). 



(2) 



In this expression, A 2 is the hyperangular kinetic energy; 
and V, the potential energy. The eigenfunctions $„ of 

H a d, 



H ad (R,n)<Z>„(R;{l) = U v {R)$ v (R;Sl), 



(3) 



form a complete, orthonormal basis at each R. Tthe total 
wave function ip{R, f2) can thus be written as 



(4) 



II. THEORETICAL METHOD 

Our method for solving the three-body Schrodinger 
equation has been detailed elsewhere 0, HH , but we in- 
clude a brief description here for completeness and clar- 
ity. 



A. Adiabatic hyperspherical representation 

We solve the three-body Schrodinger equation using 
the adiabatic hyperspherical representation [24], [2(| . Af- 
ter the usual separation of the center-of-mass motion, 



where v is a collective index that includes all the quan- 
tum numbers necessary to identify each channel. If the 
expansion in Eq. Q includes the complete, countably- 
infinite set of ip(R,Cl), then this representation of 
if) is in principle exact. In practice, of course, the sum 
is truncated to a finite number of terms, but it can be 
extended systematically to obtain essentially any desired 
level of accuracy. The eigenvalues Uy{R) of Eq. © give 
the potential curves. 

Substituting Eq. (TJ| into the full Schrodinger equa- 
tion (fT]) and projecting out $^ gives the hyperradial 
Schrodinger equations — a system of coupled ordinary 
differential equations: 



J 
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U v (R)^j F V {R) -1-J2( 2P ^'( R )^ + Qw>(R)) Fw(R) = EF V {R). 



(5) 
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The nonadiabatic coupling terms P vv i(R) and Q VV >(R) 
are generated by the hyperradial dependence of the chan- 
nel functions and are responsible for inelastic transitions. 
They arc defined as 



P UV ,(R) = (($„(#) 



d 

dR 



and 



Q^(R) = (($ V (R) 



dR 2 



(6) 



(7) 



The double brackets indicate integration only over the 
angular coordinates Q and traces over any spin degrees 
of freedom, with the hyperradius R fixed. 

The S'-matrix elements — and thus cross sections and 
rates — are found from the solutions of Eq. §5§ [13] ■ The 
main focus of this work, the three-body recombination 
rate K3, is obtained from the S'-matrix as follows fllll24l|: 



EE 

J,7T I , f 



192(2 J + 1)tt 2 
JH 4 



(8) 



where k = ^/2jlE is the incident hyperradial wave num- 
ber, i labels the initial continuum channel, and / labels 
the final two-body channels. This expression also reflects 
the fact that we use the total orbital angular representa- 
tion since the total orbital angular momentum J is con- 
served. Total parity ir is also a good quantum number 
labeling the S-matrix. 



B. Two-body potential model 

The potential V used in Eq. ([2]) for the present cal- 
culations is a sum of atom-atom interactions, V = 
v{rx2) + v(r3i) + v(r23), which is most appropriate for 
spin-stretched atoms. We take advantage of universality 
to use a two-body model potential for v(rij) that is com- 
putationally convenient. That is, the Efimov effect and 
other low-energy three-body properties do not depend 
on the details of the interatomic interaction, only on the 
scattering length and the characteristic length scale r of 
the potential. The primary simplification this allows is 
to reduce the number of two-body ro-vibrational bound 
states from hundreds or thousands, for realistic triplet 
alkali potentials, to something more manageable. 

The two-body potential model adopted here is 



v(r) = D sech 2 



(9) 



where D is the potential strength. With this potential, 
we can easily produce two-body systems with different 
values of the scattering length and different numbers of 
bound states by changing D 0, [ll[ . 

In Fig. [T] we show the scattering length as a function of 
D. The first pole in a (as D grows increasingly negative, 



starting from D = 0) occurs when the first s-wave bound 
state is formed. The second pole in a (as \D\ increases) 
occurs when a second s-wave bound state is formed. Our 
rate calculations cover regions I and II. In addition to 
having one more s-wave bound state than region I, re- 
gion II also has a <i-wave two-body state that becomes 
bound at the value of D indicated in the figure. As we 
will see below, this d-wave state leads to resonance ef- 
fects in the three-body recombination rate that could be 
misinterpreted as an Efimov feature. 
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FIG. 1: Scattering length as a function of the potential depth 
for v(r) from Eq. (0. There is one s-wave bound state in 
region I but there are two in region II. A d-wave bound state 
is also created in region II. 



III. RESULTS AND DISCUSSION 

The fact that three-body recombination can be the 
main loss mechanism in ultracold atomic gases is what 
made Efimov physics observable in the recent experi- 
ments of Kraemer et al. For a > 0, the recombination 
rate at T—0 shows a series of minima, resultin g fro m in- 
terference between two different pathways [TTL Il2|. The 
rate can be conveniently written as 



Ko = 67.1e- 2m 



sin 2 [s ln( — ) + $ M ] + sinh 2 r] M 
ro 



fin 4 

16.8(1 - e - 4 ^) — , 
m 



ha 4 
m 

(10) 



including the contributions from both weakly- and 
deeply-bound molecules @, [ll|, EH- F° r a < 0, 
though, the mechanism that produces signatures of Efi- 
mov physics is quite different and is related to resonant 
transmission effects that occur when an Efimov state is 
created [ll[. F° r a < 0, the T~0 recombination rate is @ 



4590sinh(277p) 



H\a\ 



sin [sq \n(\a\/r ) + $p] + sinh T] P 



(11) 



In these equations, $m and ^p are short-range three- 
body phases that determine the positions of the minima 
and peaks, respectively. The parameters r\M and r\p were 
introduced Q to characterize the probability of an in- 
elastic transition at small distances to a deeply bound 
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molecule. In practice, $ and 77 are used as fitting param- 
eters since the short distance behavior of realistic systems 
is not generally known. We have found Q that Eqs. (fT0|) 
and are good approximations at all collision energies 
E and scattering lengths a in the threshold regime, i.e. 
E < ft 2 /2/i2ba 2 where /X2b is the two-body reduced mass. 



A. Fitting experiment 

In Fig. [Ha) we show our numerical results for the 
thermally averaged @, [H| recombination length £3 of 
133 Cs atoms, along with the experimental data obtained 
in the Innsbruck experiment [5]. Our data were gener- 
ated for the same range of scattering lengths (in region 
I of Fig. [1]) and temperatures accessed experimentally, 
including the two leading-order contributions, J n = + 
and 2 + , to Eq. §8§ [7\. Even though we used a sim- 
ple two-body model potential, our results for a < and 
a > reproduce the experimental data quite well. In 
our approach, the only adjustable parameter is tq. We 
obtained ro~100 a.u. by fitting the experimental peak 
position, a w —850 a.u., at T—10 nK. We used this same 
ro to produce the curves for the other temperatures and 
for a > 0. We will discuss the validity of the latter in 

Sec. Unci 

It is important to note that although our numerical 
value for ro happens to agree with the expected result for 
Cs atoms — namely, the van der Waal's length r v nw 
— it unfortunately does not imply that we can predict the 
peak position for this or any other species. The reason 
is that our choice of model two-body potential Eq. © 
determines the short-range physics and has no particular 
relation to the correct short-range physics for three Cs 
atoms. A different choice of model potential would lead 
to a different best fit value for r$. The model potential 
thus sets the values of $ and 77, preventing their use as 
fit parameters. In fact, any finite range two-body model 
potential [20|, |2l| will suffer from essentially the same 
problem, as we will discuss in more detail in Sec. IIIID1 

The major difference between our results in Fig. [Ha) 
and the experimental results lies at a > 0. As it turns 
out, we would not expect good agreement since our a > 
results use the same r as the a < results. As already 
noted, this point will be discussed in detail in Sec. IIII CI 
but the a > results in Fig. [5fa) do serve to illustrate 
the important role of temperature for attempts to ob- 
serve the minimum. Figure[^a) shows, for instance, that 
at 10 nK we do find a minimum at a — 1500 a.u.. At 
250 nK, however, strong contributions from the next par- 
tial wave, J n = 2 + 7], mask the minimum. Comparison 
with the experimental data shows that even if there were 
actually a minimum for Cs that matched our model, the 
experimental temperature was likely just too high to be 
able to see it. For heavy atoms like Cs, the requirement 
of being in the threshold regime to see Efimov features [7[ 
places a relatively severe limit on the maximum allowable 
temperature. 



B. Finite range potentials 

Analytic expressions for K3 such as Eqs. (fTU|) and 
(fTTj) were derived as an expansion about \a\ ^> r . In 
practice, this was accomplished using zero-range two- 
body interactions that support at most one s-wave bound 
state. Calculations performed using a finite range poten- 
tial like Eq. (J9]) produce deviations from these analytical 
predictions and introduces features not observable with 
these zero-range potentials. One immediate consequence 
of a non-zero ro is that the minimum observed in the 
Cs experiment around a = 210 a.u. is likely not re- 
lated to Efimov physics since a is not much larger than 
= TvDW ~100 a.u. This point was, in fact, already 
noted in Ref. [j| . We have found § that for most quanti- 
ties the predictions of the zero-range model are recovered 
in a reasonably quantitative way only when \a\ is more 
than an order of magnitude larger than tq. 

In this range of scattering lengths, the recombination 
rate for finite range potentials can have structure that 
is completely independent of Efimov physics [l7| . This 
structure is generated by the presence of non-s-wave two- 
body states. Figure [2jb) shows an example of such struc- 
ture due to a d-wave two-body state. It shows our nu- 
merical results for Q3, now obtained for a in region II 
(with two s-wave bound states). These results display 
behavior similar to the experimental observations, rep- 
resented in this figure by the loss fraction reported in 
Ref. Q . A minimum (or peak) in the recombination rate 
does lead to a minimum (or peak) in the loss fraction, 
and vice-versa, so the figures can be compared quali- 
tatively. As the scattering length is decreased from its 
maximum in the figure, the two-body d-wave state can 
be visualized as approaching the threshold from above. 
While above the threshold, the d-wave state is a two- 
body shape resonance; below threshold, it is, of course, 
a bound state. Thus, the first, sharp feature as a is de- 
creased, at a w 280 a.u., is produced when the two-body 
shape resonance energy matches the three-body collision 
energy [TJ]. Physically, this situation can be thought 
of as ro-vibrational relaxation of the two-body resonant 
state. Two atoms collide, forming a d-wave two-body 
resonance, when the third atom collides with it, driv- 
ing it down into the available s-wave bound states. It is 
interesting to note that the total orbital angular momen- 
tum is J = 0, whereby the relative angular momentum 
of the initial atom-dimer complex must also be a d-wave. 
The final angular momentum of the bound dimer is zero 
as must be the relative atom-dimer angular momentum. 
There is thus considerable angular momentum and en- 
ergy exchange at work near this feature. The next broad 
feature at smaller a occurs where this two-body resonance 
becomes a bound state. Being the highest-lying bound 
state, it receives the majority of the recombination. 

It is worth noting that this d-wave K3 resonance does 
not correspond to a resonance in a because it is an open 
channel two-body shape resonance — not a resonance in 
the s-wave channel. Higher partial wave resonances in a 
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FIG. 2: (Color online) (a) Comparison between 133 Cs recombination length, £3 = (mK3/h\/3) 1 ^ 4 ' [ill ], obtained experimentally 
[H (black symbols) and our present thermally averaged results (red symbols). Filled circles and triangles were obtained at 10 
and 200 nK, respectively; open diamonds and squares were obtained at 250 nK and 250-450 nK, experimentally, and at 250 nK 
in the present calculation, (b) Effect of the d-wave resonance on the thermally averaged recombination length for small values 
of a in region II of Fig. [1] in comparison to the observed loss fraction from [5j. 



are typically due to higher order two-body interactions 
couplin g t he incident s-wave to a closed channel reso- 
nances pjf . The appearance of higher partial wave reso- 
nances at moderate scattering lengths is a generic feature 
that has been examined in detail for van der Waals inter- 
actions in Ref. 16]. The two-body interaction for most 
any atoms will generally support not only d-wave bound 
states and resonances, but also much higher angular mo- 
mentum states, potentially leading to considerably more 
structure near a = 0. 

We emphasize that the comparison in Fig. HJb) is 
only an illustration of the complications that can oc- 
cur in realistic systems, especially in the non-universal 
region of \a\ « tq. A definitive interpretation of the 
a > experimental data requires incorporating consider- 
ably more information about the Cs-Cs interaction into 
our model than we have. Fortunately, much is known 
about the Cs-Cs interaction, and the scattering length 
has been calculated as a function of magnetic field with 
some confidence [Hj]. With that information, and the 
additional universal relations between s-wave and e?-wave 
states for alkali atoms [lj|, it should be possible to lo- 
cate any higher partial wave open-channel resonances to 
check whether features like the d-wave K3 resonance in 
Fig. [2fb) are possible. 

Further complicating any effort to identify Efimov fea- 
tures in this range of a, the global minimum in K3 lies 
at a ~ 200 a.u. — even without the d-wave resonance. 
The global minimum can be see in both Figs. Hfb) and 
(c) , showing that recombination goes through a minimum 
when a is tuned from a — > 00 and a — > — 00 through 
a = 0. 

To further illustrate the effect of finite range two-body 
potentials, we show in Fig. [3] the recombination rate ob- 
tained for fixed collision energies from 0.67 x 10~ 4 nK up 
to 670 nK and \a\ up to 2 x 10 6 a.u. across regions I and 



II of Fig. [T] These rates include only the + contribu- 
tion and are not thermally averaged, so the indicated 
"temperatures" are the collision energies converted to 
Kelvin. The rates in Figs. Ela) and [3fb) were obtained 
for a in region II and therefore have up to three con- 
tributions: recombination into the two possible s-wave 
states and into the <i-wave state. The rates in Figs. \>2,c) 
and[3]Jd) were obtained for a in region I and represent re- 
combination into the single available s-wave bound state. 
Figures [3^a) and (b) thus connect across a zero in the 
scattering length; (b) and (c), across a pole; and (c) 
and (d), across a zero. Comparing the partial rates in 
Figs. [3]^b) and (c) for energies from 6.7xl0~ 2 nK up 
to 670 nK, recombination into the lowest s-wave bound 
state is smooth across the resonance where the second 
s-wave bound state is formed. For lower energies this 
connection should also be smooth, but will happen only 
for larger values of a. Recombination into the second s- 
wave bound state in Fig. \Sj[h) introduces a discontinuity 
in the total recombination rate at non-zero collision ener- 
gies. It would apprear from Fig.[3Jb) that the appearance 
of the ci- wave bound state near a=0 also introduces a dis- 
continuity, but because the relative angular momentum 
between the atom and dimer in the final state is cZ-wave, 
the Wigner threshold law guarantees that it will turn on 
smoothly. 

Figure [3Ja) shows that recombination into the d-wave 
state actually dominates over the whole range of a in this 
panel, which is not too surprising since it is the most 
weakly bound dimer state. Moreover, the d-wave partial 
rate shows the expected Efimov features — as do all of 
the partial rates in Fig. OHa) — emphasizing that the 
Efimov physics is determined for a < by the initial 
three-body state and not the final state 0, EH, E3 ■ This 
dependence is one way to understand the universality of 
K3 for a < since the initial channel is universal while 



6 



the final, deeply bound states are not. 

One other issue evident from comparing the analytical 
curves in Fig.[3Jto the calculated ones is that the expected 
spacing between Efimov features does not hold between 
the two features nearest a=0. In Figs. [3jc) and EJd), 
the spacing between the first two minima and first two 
peaks are 28.3 and 14.3, respectively, instead of the 22.7 
expected for identical bosons. The spacing between the 
next two minima and peaks are 23.4 and 21.6, respec- 
tively. Where the spacing can be defined in Figs. [3ja) 
and (b), it follows the same pattern. This sort of devia- 
tion for the spacing of the first few features is important 
to be recognized by any experiment seeking to verify the 
logarithmic spacing of Efimov features. At least in the 
near future, experiments will likely not be able to ob- 
serve more than a few of these features nearest to a = 0, 
in view of the large range in a and the low temperatures 
required Q. So, knowing the finite range corrections to 
the logarithmic spacing, and knowing whether those cor- 
rections are universal, are both critical. 



C. Tuning a from oo to +oo 

In addition to the complications discussed above, fun- 
damental issues were raised in Ref. @ since the scattering 
length was tuned from -co to +oo through a = and 
therefore through a non-universal region. In fact, a was 
tuned through a series of narrow resonances. While this 
was obviously the expedient choice experimentally, it is 
not obvious that the a < and a > Efimov features 
should retain their universal relation predicted from the 
zero-range model [9(, 

A$ = $ M - $ P = -1.53(3), (12) 

or even that the short-range three-body phase on each 
side of a narrow resonance should be universally related, 
let alone the same. Clearly, it is convenient if these rela- 
tions hold since only one phase is then required to predict 
the position of both minima (a > 0) and peaks (a < 0). 
So, the question that naturally arises is whether there is 
a relation similar to Eq. (|12p between the peak and min- 
imum positions when a is tuned across a non-universal 
region (a = 0). To explore this issue, we analyze the re- 
lations between the numerical calculations shown in the 
panels of Fig. [31 The fact that K$ for a < and a > 
are universally related only across a pole in a, and the 
resulting consequences for the Cs experiment, have been 
raised in previous work [TcL l20j| . Here, however, we show 
via a counter-example that there is no general universal 
connection across a zero in a. 

Because of the slow approach of the system's properties 
to the \a\ 3> ro predictions and the requirement of being 
in the threshold regime, the range of scattering lengths 
and temperatures required to observe more than one fea- 
ture in 133 Cs, necessary to determine $ precisely from 
Eqs. (|10p and (fTTj) . greatly exceeds current experimental 



capabilities. This fact underscores the advantages of us- 
ing different atomic species to make the required range 
of a and T much more reasonable @ ■ 

The three-body phases <&p and $m were determined 
from our numerical results in Fig. [3J by independently 
fitting Eqs. (|10[) and (fTTj) to the positions of the peaks 
and minima at our largest \a\ and lowest temperature in 
order to better satisfy the approximations under which 
those equations were derived. The fit curves are indicated 
in Fig. [3J by the line without symbols. From Figs. [3](b) 
and (He), which should be connected universally since 
they cross a pole in a, we obtain A<&b c =— 1.553, in good 
agreement with Eq. (fT2"|) . The other phase differences — 
which imply crossing a non-universal region — are, how- 
ever, A$ ab =-1.053 and A$ cd =-1.381. Thus, the Efi- 
mov physics for a>0 and o<0 do not appear to be simply 
connected when separated by a non-universal region. 

We rationalize this non-universal behavior by recog- 
nizing that small changes in the interactions should not 
substantially affect the short range physics, i.e. $ and rj. 
Changing a from positive to negative across a resonance 
requires a much smaller change in the interactions than 
does crossing a — 0, so the former should lead to univer- 
sal behavior. The relative changes required in the depth 
of the potential D when crossing each region are clearly 
illustrated in Fig. [T] Similarly, when controlling a via a 
Feshbach resonance, going from a > to a < requires 
smaller changes in the magnetic field when crossing a res- 
onance than when crossing a — 0. It is not entirely clear, 
however, whether this argument generalizes directly from 
our model potentials to realistic potentials with many 
bound states since the relative change in such a realistic 
potential when crossing a—0 will also be small. In princi- 
ple, we can answer this question by simply continuing to 
increase D and solving the problem numerically. Unfor- 
tunately, we are not yet able to calculate the three-body 
recombination in a system with more than a handful of 
two-body bound states. Incidentally, we expect similar 
behavior for the other short-range three-body parame- 
ter r\. Specifically, we expect rjp = r\M when a is tuned 
across a resonance, but not when it is tuned across a = 0. 

In the Innsbruck experiment, A<I> was found in the 
form of the ratio 

^=cxp[-(A<f + 7r/2)/ So ]. (13) 
\a-\ 

Experimentally, a+ located the first recombination max- 
imum near a — for a > 0; and a_, the first peak for 
a < 0. The experimental fit gave o+/|a_ I w 1.25(9) H. 
Considering that the phases in A<f> are defined as mod 
7T, with $m and $p defined in the range [0, 7r), the min- 
imum and maximum values for a + /a_ are respectively 
0.00925 and 4.76385, this result agrees relatively well 
with the theoretical value of 0.96(3). Theoretically, how- 
ever, A$ is defined from the \a\ 3> ro limit. Our A<Iv, 
for instance, was determined by fitting in this limit and 
gives a+/|a_|=0.982. The proper comparison with ex- 
periment, using the first features in Figs. [3Jb) and[3]Jc), 
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FIG. 3: (Color online) Recombination rate for a covering regimes with different numbers of two-body bound states, as indicated 
above. In (a) and (b) recombination has contributions from the deeply (red circles) and weakly (green triangles) bound s-wave 
states as well as from the d-wave state (blue diamonds). The lines without symbols are the fit of Eqs. (1 10 p and (|11[) . 



is difficult to make, though, since there is no clear max- 
imum for small positive a — save for the d-wave reso- 
nance — to define a+. If there were, it would likely be 
smaller than predicted by Eq. (TIT)]) (compare with the 
solid line fit). On the other hand, a_ would be larger 
than predicted, giving o+/|a_| smaller than the \a\ ^> r 
prediction — assuming these feature shifts are universal. 
We can also get a sense of the quality of the agreement 
between theory and experiment by calculating a+/|a_| 
from A$ across a non-universal region: A$ c ^, for in- 
stance, gives a + /|a_|=0.828; A<I> ab gives 0.598. Both of 
these values are nearly as close to the predicted value as 
the experiment, but we know that the present "agree- 
ment" is accidental. So, the fact that the experiment 
accessed the non-universal region together with the fact 
that the minimum position (a = 210 a.u.) is not firmly in 
the universal regime (a ^> ro), make us believe that the 
agreement between experiment and theory for a+/|a_| is 
likely fortuitous. 



D. Three-body short-range physics 

It should be clear by now that knowledge of the short- 
range three-body phase is essential to make quantitative 
theoretical predictions for ultracold three-body recombi- 
nation of realistic systems — or, indeed, for nearly any 
ultracold three-body process. The three-body phase can 
be determined from essentially any ultracold three-body 
observable, and many examples have been discussed in 
Ref. [9j. It should be equally clear from the discussion 
of Fig. [3| that fitting a two-body model potential to give 
the correct scattering length is not sufficient to determine 
the three-body phase. Even a more c omp lete fitting of 
the low-energy two-body physics [20|, |2lJ will not give 
the correct three-body phase. In fact, quantitative pre- 



model calculations such as 



dictions for real systems would not be possible even if 
the complete two-body interaction were included in the 
calculation. The reason is simple: in real triatomic sys- 
tems, there is a short-range, non-additive, purely three- 
body interaction [22|, |23|. S o 
ours and those in Refs. [20|, [2l| , 

that fit only two-body physics fundamentally can- 
not predict the positions of Efimov features quantita- 
tively without the input of a known near-threshold three- 
body observable — because of the three-body short-range 
physics. 

To illustrate the effects of non-additive three-body in- 
teractions, we have performed calculations that include 
one such correction. We have included the well-known 
three-body Axilrod- Teller potential [22J which, analogous 
to the two-body van der Waals interaction, is the leading- 
order dispersion term describing the long- range in- 
duced dipole-dipole-dipole interaction among three neu- 
tral atoms. The Axilrod- Teller potential assumes the fol- 
lowing form, 



WAT{ri2,r23,ri 3 ) = 7 



(1 + 3 cos 812 cos 823 cos 831 ) 



,(14) 



12' 23' 31 



where are the interatomic distances and 0y = 
cos (r^ •fjk) are the inner angles of the triangle formed 
by the three atoms. In the above equation 7 is a posi- 
tive quantity whose exact value depends on the details 
of the three-body system. Note that this interaction can 
be attractive or repulsive depending upon the shape of 
the three atom system. It is also worth emphasizing that 
the Axilrod- Teller interaction is just one of several short- 
range three-body interactions that can affect three-body 
observables. In fact, much stronger corrections are ex- 
pected due to purely three-atom exchange effects [23| . 

While 7 in Eq. (] 14f) normally reflects the details of the 
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a(1000a.u.) a(1000a.u.) 

FIG. 4: (color online) Comparison of the (a) peak and 
(b) minimum position in the recombination length with the 
Axilrod-Teller interaction (green filled circles) and without it 
(red open circles). These calculations demonstrate that the 
main effect of the inclusion of the Axilrod-Teller potential is to 
move the position of the minimum and peak as shown above. 

atoms, for the purposes of our model, we have chosen it 
to have the form 

7 = 7ol-Dkn[tanh(— ) tanhf— ) tanh(^i)] 3 . (15) 
r o r r 

This choice ensures that Wat scales as we tune the two- 
body scattering length, that it scales properly with ro, 
and that the r^ 3 singularities are cut off near = 0. 
We have arbitrarily chosen a value for 70 that causes a 
variation of about 15% in the adiabatic hyperspherical 
potentials near the minimum. This magnitude change is 
probably an underestimate, given that it has been shown 
that non-additive three-body terms make the minimum 
in the three-atom potential surface from 60% to a factor 
of four deeper for alkali atom systems, compared to the 
pairwise additive potential surface [3lT I32I |33| . 

Figure [¥] shows the recombination length with and 
without the Axilrod-Teller term. Our results in the fig- 
ure demonstrate that the minimum and peak positions 
are clearly shifted by the inclusion of the Axilrod-Teller 
term. For a < 0, besides changing the peak position, 
the Axilrod-Teller interaction also changes the width of 
the resonance peak and the overall amplitude of recom- 
bination length. This result indicates that the coupling 
to deeply bound decay channels is also affected by the 



Axilrod-Teller interaction. The analytic expressions in 
Eqs. (HHJ) and (fTTj) account for these effects through 
changes to the short-range parameters $ and 77. 

Our results in Fig. [H therefore, demonstrate our ini- 
tial statement. The inclusion of purely three-body inter- 
actions such as the Axilrod-Teller — or the three atom 
exchange — are of crucial importance for accurately pre- 
dicting features related to Efimov physics. As a conse- 
quence, even though previous works have rather carefully 
fit the two-body physics [13, Hl[ , it is unlikely that they 
have predictive power in locating Efimov features. 
IV. SUMMARY 

Using essentially exact numerical solutions for a model 
two-body potential, we have demonstrated the need for 
caution when applying existing analytic expressions for 
three-boson recombination to experiment. While they do 
give considerable insight, their predictions for the first 
features near a=0 are not quantitative and may even fail 
qualitatively due to the complications arising from real- 
istic potentials — in particular, from their finite range. 
Moreover, care must be taken to ensure that the tem- 
perature is truly in the threshold regime. We have also 
verified quantitatively that a > features are universally 
connected to a < features only when a is tuned through 
\a\ = 00 of a given resonance. 

Finally, we have emphasized the importance of the 
short-range three-body phase in making quantitative pre- 
dictions for real systems. This point does not appear to 
be fully appreciated theoretically, but our example of the 
non-additive three-body term underscores the limitations 
of purely two-body models for actual triatomic systems 
and thus the need for a three-body parameter. Unfor- 
tunately, purely ab initio calculation of this three-body 
parameter seems, for now, out of reach for the systems 
of experimental interest. 
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